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This document was submitted as supporting material to an Engineering Change Proposal (ECP) 
for the Square Kilometre Array (SKA). This ECP requests gridded visibilities as an extra imaging 
data product from the SKA, in order to enable bespoke analysis techniques to measure source 
morphologies to the accuracy necessary for precision cosmology with radio weak lensing. We also 
discuss the properties of an SKA weak lensing data set and potential overlaps with other cosmology 
science goals. 


I. RADIO WEAK LENSING COSMOLOGY 

Probing cosmology with the weak gravitational lensing effect is potentially the most sensitive way of determining 
the evolution of the Dark Energy equation of state w with redshift and hence determining its physical nature. Weak 
lensing is therefore a primary science driver for high-profile upcoming optical experiments such as Euclid and LSST. 
However, weak lensing requires exquisitely precise measurements and a number of irreducible systematic effects may 
ultimately limit the power of these optical surveys. Large, deep surveys with SKA-MID will also be capable of 
achieving the number density of high redshift sources necessary for weak lensing, with both SKAl and the full SKA 
(SKA2) providing constraints on cosmological parameters competitive with optical experiments, as shown in Fig. 

In addition, there are a number of unique advantages to performing weak lensing in the radio (see [3] for more detail): 

• Highly stable and deterministic point spread functions. 

• Source redshift distributions with higher median redshifts, giving more constraining power on evolution of Dark 
Energy and other cosmological parameters over cosmic time. 

• Extra information from polarisation and spectral line velocity measurements which may be used to remove 
otherwise limiting astrophysical systematics such as intrinsic galaxy alignments. 
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Figure 1. Forecast constraints on total matter content Hm and clustering as from stage III (left) and stage IV (right) weak 
lensing experiments, assuming source shapes can be measured to the same accuracy in different wavebands. 
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• Cross correlations between optical and radio leasing maps may be expected to have uncorrelated telescope 
systematics, allowing the best unbiased constraints. 

However, weak leasing has stringent and specific requirements on image fidelity; the cosmological shear signal which 
affects source ellipticity is small and so must be averaged over many source shapes. This means source ellipticities 
must be measured with a systematic error of less than 1 part in 10,000 in order for errors on cosmological parameters 
to be dominated by statistical, rather than systematic, errors. 

Simulations have shown that, even on high SNR objects (whose shapes are easiest to measure) with an SKAl 
antenna configuration, current iterative deconvolution methods produce images with structures in the residuals which 
dominate the cosmological signal, as shown in Fig. from [5]. The forecasted constraints displayed in Fig. [^assume 
that shapes can be measured with equivalent accuracy in the radio as is currently possible in the optical. However, 
using ellipticity measurements with the properties of those in Fig. would produce constraints with analysis-induced 
biases more than 100% of the target measurement of ic = — 1, far larger than the achievable 1% statistical uncertainty. 
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Figure 2. Biases on ellipticity from CLEAN images of a simulated SKAl-MID data set from [5]. Top: input ellipticity plotted 
against difference to recovered ellipticity, showing the considerable multiplicative (slope) and additive (offset) biases. Bottom: 
requirements on biases for a 5000 deg^ SKAl survey to have errors dominated by statistics, along with values recovered from 
CLEAN images. 

We stress that this bias is not a fundamental limitation of the SKA telescope, which will be capable of capturing all 
of the information necessary for precision weak leasing, but is induced by the analysis procedure adopted for turning 
the visibility data into images. This ECP requests that gridded visibilities be made available, which will facilitate 
analyses more optimal for the particular needs of weak leasing, thus enabling this frontier science. 

The current SKA Phase 1 Level 1 requirements document provides for: moment images of a multi-frequency 
synthesis imaging process as the data products from the continuum imaging pipeline; and spectral line and continuum 
model images from the spectral line emission pipeline. In this ECP we request as an additionaly available data 
product the gridded visibilities which constitute the Fourier-plane representation of the dirty map, prior to the final 
imaging step. This will allow use of both bespoke imaging techniques for weak lensing, which may be tailored to the 
particular needs of this science case, and the potential for directly modelling the visibility plane data (as was done in 
the detection of a radio weak lensing signal by 0 ). 


II. REQUIREMENTS FOR SKA 

Whilst gridded visibilities do not represent the full amount of information available from the raw visibilities, the 
operations performed when gridding should be either directly computed as part of the analysis pipeline (e.g. sampling, 
weighting, convolution with gridding kernel) or precisely known from features of the telescope (convolution with A and 
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Figure 3. Requirements on bandwidth (left) and time (right) smearing for systematic errors on measured ellipticity to be below 
statistical errors for SKAl and SKA2 weak leasing surveys at the edige of the primary field of view (25 arcmin from the pointing 
centre). 


W-projection kernels), meaning their effect on source ellipticities should be able to be modelled. This is in contrast 
to the effect of an iterative deconvolution imaging algorithm on an unknown sky. 

Requirements on these steps which move from raw to gridded visibilities will take the form of adequate samplings 
of the spatial scales of interest for weak lensing source modelling; adequate removal of bright confusing sources in 
far sidelobes; adequate w-projection and adequate primary beam calibration; and ensuring relevent morphological 
information is not destroyed by time and bandwidth smearing (as discussed below). Whilst it may be the case that 
many of these requirements will be satisfied along with other science cases requiring high dynamic range imaging, we 
expect to carry out detailed simulations to confirm exact requirements on the levels of sampling required. However, 
even within current default specifications for the gridding process, the ability to perform tailored imaging and source 
modelling analyses from gridded visbilities will be essential for enabling precision cosmology with radio weak lensing. 

In order to get a sense of the necessary data set size for a weak lensing survey, we calculate the effect of time and 
bandwidth smearing on the ellipticities of sources across the field of view. We follow [2] and calculate the change in 
the FWHM of the sythesised beam in the radial (tangential) direction for bandwidth (time) smearing at the edge 
of the field of view (25 arcmin from the pointing centre) for an SKAl-MID experiment with a circular synthesised 
beam FWHM at the pointing centre of 0.5 arcsec. We vary the channel width and compare the maximum induced 
ellipticity (at the edge of the field of view) to requirements on the systematic uncertainty (as described in [3]) on PSF 
ellipticity for SKAl and SKA2. The result is shown in the left panel of Fig.[^ finding ^ 50 kHz channels are required 
for the bandwidth smearing to be tolerable; for an observation making use of 30% of the total available bandwidth 
(i.e. 300 MHz), this means ^ 6000 channels will be required. Similarly for the time smearing, the right panel of Fig.[^ 
shows integration times ^ 1 second will be required for smearing-induced ellipticity to be at acceptable levels. 

Whilst the smearing from both of these effects should be deterministic and hence correctable, assessing the residuals 
from such corrections across the full space of observations would require full-scale (i.e. on the same scale as the data) 
visibility plane simulations. In this case, the requirements on smearing would then be in terms of required residuals 
after correction (which will likely be limited by the variations in smearing across the size of a typical source). However, 
generating and processing the massive simulations which such a correction scheme would require represents a huge 
computational challenge, equivalent to the collection and processing of the real survey data. It therefore makes much 
more sense to place requirements on the averaging procedures applied during the data processing, thus ensuring that 
these systematic effects are suppressed to acceptable levels in the initial data products without further corrections. 


III. SUMMARY 

We expect this request to necessarily require no new algorithm development on the part of the SKA Organisation 
alone, with this being the responsibility of the science case. The impact on the telescope operations are expected to 
be as follows: 











4 


• An addendum to the spectral line emission or continuum data pipeline which outputs as a data product the full 
gridded visibility cube (before FFT to a dirty image cube and deconvolution) 

• Requirement for a larger data archive for observations comprising the weak lensing survey. For a 5000 deg^ 
survey to a depth of 2 /iJy RMS, we expect to use ^ 10^ pointings of ^ 1 hours each, requiring ^ 6000 frequency 
channels at a resolution of 50 kHz. This is a large multiplicative facor in data volume compared to a continuum 
survey, but is directly comparable to that expected by HI line galaxy surveys, who specify 10 kHz channels over 
similar total bandwidth and total observing time [5] . 

• Alternatively^ a weak lensing dedicated pipeline, or sub-pipeline, which is capable of measuring source shapes in 
cached visibilties to the required level of accuracy. 

• A corresponding increase in data transport from the Science Data Processor to the science analysis team. 

• An increase in complexity for the Telescope Management system to allow for this mode of data output. 

We expect this weak lensing survey to have similar requirements in terms of time, depth and footprint to other large, 
deep surveys at L-band (such as SKAl Science Goal 27: “The resolved all-Sky characterisation of the interstellar and 
intergalactic magnetic fields”, SKAl Science Goal 33: “Angular correlation functions to probe non-Gaussianity and 
the matter dipole” and SKAl Science Goal 36: “Test dark energy & general relativity with fore-runner of the ‘billion 
galaxy’ survey.” [T]), other users of which will benefit from the availability of gridded visibilities as a data product. 

The SKA telescope will have the capacity to revolutionise the field of weak lensing by extending the reach of this 
key cosmology experiment to the radio band. Not only will the SKA be competitive with state-of-the-art optical 
experiments (such as Euclid and LSST) in terms of raw precision, it will offer truly unique approaches to weak 
lensing (that are not possible with Euclid or LSST) through polarisation and HI rotational velocity measurements. In 
order to realise the huge science return available, bespoke galaxy shape analyses will need to be applied to the data. 
The gridded visibilities and/or weak lensing optimised sub-pipeline that we request in this ECP will facilitate the 
application of such bespoke analyses, thus ensuring that precision weak lensing cosmology remains within the grasp 
of the SKA. 
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